Nonuniversal finite-size scaling in anisotropic systems 
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We study the bulk and finite-size critical beliavior of tlie 0(n) symmetric theory with spatially 
anisotropic interactions of non-cubic symmetry in d < 4 dimensions. In such systems of a given (d, n) 
universality class, two-scale factor universality is absent in bulk correlation functions, and finite-size 
scaling functions including the Privman-Fisher scaling form of the free energy, the Binder cumulant 
ratio and the Casimir amplitude are shown to be nonuniversal. In particular it is shown that, for 
anisotropic confined systems, isotropy cannot be restored by an anisotropic scale transformation. 
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A basic tenet in the physics of critical phenomena is 
the notion of a universality class. It is characterized by 
the dimensionality d of the system and by the number n 
of the components of the order parameter. (See, e.g., the 
review article flj.) Within a certain {d,n) universality 
class, the universal quantities (critical exponents, am- 
plitude ratios and scaling functions) are independent of 
microscopic details, such as the particular type of (finite- 
range or van der Waals type) interactions or the lattice 
structure Q ■ This implies that a given universality class 
includes both spatially isotropic and anisotropic systems. 

Once the universal quantities of a universality class 
are known the asymptotic critical behavior of very dif- 
ferent systems (e.g., fluids and magnets) is beheved to 
be known completely provided that only two nonuniver- 
sal amplitudes Ai and A2 are given. This property is 
known as two-scale factor universality or hyperuniversal- 
Q> 0> 0| • In terms of the singular part of the reduced 
bulk free energy density Fs/VksT = fs{t, h) above Tc, 



f,{t,h) = A^t^-W{A2ht~P') 



(1) 



with PF(0) = 1 and t = (T - Tc)/Tc < 1, this property 
can be stated as 

^lirn^ fs{t-, 0)C'' = Q{d, n) = universal . (2) 

Thus the amplitude ^0 ~ (Q/^i)^^'' of the correlation 
length ^ = Coi at zero ordering field h is not an inde- 
pendent amplitude but is universally related to Ai . The 
validity of two-scale factor universality has been estab- 
lished by the renormalization-group (RG) theory on the 
basis of an isotropic Hamiltonian with short-range inter- 
actions below the upper critical dimension d* = 4 01 but 
no general proof has been given for the anisotropic case. 

In this paper we study the critical behavior of sys- 
tems with a spatial anisotropy of non-cubic symmetry 
within a given (d, n) universality class. An example is 
an Ising ferromagnet with an isotropic nearest-neighbor 
(NN) coupling J > and an anisotropic next-nearest- 
neighbor (ANNN) coupling J' on a simple-cubic lattice. 
In some range of J' / J this model has the same type of 



critical behavior as the ordinary (J' = 0) Ising model. 
We shall show that for such systems Eq. ^ must be 
generalized to 

d 

lim /a(i,0) = Q{d,n) = universal (3) 

1=1 

where — Co*''^"'' are the correlation lengths associated 
with the principal directions of the anisotropic system 
and where Q{d, n) is the same universal quantity for both 
isotropic and anisotropic systems. (For d = 2,n = 1, this 
is already known for the Ising model with anisotropic 
NN (ANN) interactions 7^ Jy 3].) There are, in gen- 
eral, d+1 nonuniversal bulk amplitudes ^q^"* , . . . , ^q'*'* , A2 
whose ratios are also nonuniversal. Note that there still 
exists a unique critical exponent i'{d, n) that is identi- 
cal for isotropic and anisotropic systems within the same 
(d, n) universality class 6, 7, 8, 9, lOj. 

A different type of critical behavior exists in the so- 
called strongly anisotropic systems 0, lisl where 
not only amplitudes depend on the spatial directions but 
also the critical exponents (e.g., i^n and vj^) depend on 
the direction. These systems do not belong to the (d, n) 
universality class of ordinary critical points and our anal- 
ysis will not include such types of anisotropy. 

While Eq. Q is a natural generalization of Eq. Q 
we shall call attention to the intriguing problem of finite- 
size effects in anisotropic systems. For simplicity we shall 
confine ourselves to the case of periodic boundary condi- 
tions in rectangular Li x L2 x ... x Ld block geometries 
(including L'^ cubic geometry and oo''"^ x L film geom- 
etry). There have been several studies of this problem 
in the past H E El H [H [13 . We shall only briefly 
comment on the more complicated case of anisotropic 
confined systems with non-periodic boundary conditions 

iLilEflEllillil 

It has been hypothesized p that two-scale factor uni- 
versality holds not only for bulk systems but also for con- 
fined systems, except that the finite-size scaling functions 
depend on the geometry and on the boundary conditions. 
For example, for a system in a cube of volume L"^ with 



2 



periodic boundary conditions, the singular part of the 
reduced free energy density fs{t,h,L) near bulk was 
predicted to have the asymptotic scaling form for large 

Mt, h, L) = L-'^Y.^UCitL^''' . C2hL'''/n (4) 

where the function Ycube{x,y) is universal and where Ci 
and C2 are the only nonuniversal parameters. A simi- 
lar ansatz was made for the correlation length in a 
L'^"^ X 00 cylinder j^. As a consequence, the aniplitude 
i^cuhe(0, 0) and the Binder cumulant ratio jll 124 I2,'tI| 



(at h = 0) 



U 



(d^Y^ube/dy^) I {d'Y,ube/dy^ 



y=0,x = 



(5) 



are predicted to be universal. (For example, they should 
be independent of the ratio J' /J.) The scaling form 
if extended to realistic geometries and boundary 
conditions, has far-reaching consequences for measurable 
quantities P, |2^ . In particular the prediction of a uni- 
versal character of the Casimir amplitude 



A={d-l)YfUrn{0,0) 



(6) 



is of interest, e.g., for fluid j23, superfluid |23|, and su- 
perconducting 29] films. 

The universality of the scaling functions Y of Eqs. Q 
or lO was supposed to be valid for all systems in a given 
universality class 0, including anisotropic lattice sys- 
tems provided that an appropriate rescaling of the lattice 
spacings (or length L) is performed j^. This appears to 
be consistent with existing studies of finite-size effects 
in anisotropic systems where it was stated that isotropy 
can be restored as ymp totically b y a n anisotropic scale 
transformation [ll [ll [ll jll il ii 113 . 

We have found that this picture of finite-size effects 
in anisotropic systems, though valid in special cases, is, 
in general, not correct. In the present paper we show, 
for periodic boundary conditions, that Eqs. Q - JH)), 
though valid for isotropic systems and for syst ems with 
cubic symmetry in the range i/^ < 0(1) ISlJ, are not 
universally valid for the anisotropic systems of the type 
described above (e.g. spin models with NN and ANNN 
interactions on simple-cubic lattices) although these sys- 
tems belong to the same universality class as isotropic 
systems. In such anisotropic systems the finite-size scal- 
ing functions depend, in general, on additional nonuni- 
versal parameters (apart from Ci and C2), even after a 
rescaling of the lattice spacing or of the length L. Thus, 
in general, two-scale factor universality and isotropy can- 
not be restored and the notion of a universality class is 
only of restricted relevance for the scaling functions of 
confined systems. 

We shall prove our claims within the 0(n) symmetric 
ip^ field theory with the spatially anisotropic Hamiltonian 



H{ro,uo,A;A;V;ip) = 



a, 13 



Agf} dip dip 
2 dxa dxp 



2^2 



uo('^ ) 



(7) 



for the n-component field <^(x). The sum runs over the 
components Xa of the spatial coordinates x, a = 1 , . . . , d. 
The dxd anisotropy matrix A = (A^^) is assumed to be 
real, symmetric and positive definite. This model has a 
critical point at some value rg = r'oc(A; uq, A) where A is 
a (sharp or smooth) cutoff in k space. In addition to the 
three parameters tq , uq and A of the standard isotropic 
(A = 1) model, our model has d{d+ l)/2 nonuniversal 
parameters contained in the matrix A. Below we shall 
argue that the non-diagonality of the anisotropy matrix 
A is a generic case of real anisotropic systems. For sim- 
plicity we assume a cubic volume, V ^ L'^ , Q < Xq, < L, 
with periodic boundary conditions. 

First we prove that the model defined by Eq. be- 
longs to the same bulk universality class as the standard 
isotropic Landau- Ginzburg- Wilson model with A = 1. 
The characteristic properties of the matrix A are de- 
scribed in terms of the d eigenvalues > and eigenvec- 
tors defined by Ae^ = A^e^. A rotation by the orthog- 
onal matrix U yields the diagonal matrix UAU^^ = A 
with diagonal elements Ai. After the transformation of 
the spatial coordinates 



= A 



-1/2 



Ux 



(8) 



and of the field 

(^'(x') = (det A)i/V(U"' A^/V), (9) 



det A = J]^ A, > , 

the Hamiltonian {Tj) becomes 

i7(ro, uo, A; A; V; ^) = H'{ro, u'^, A'; V; ip') = 



d'^x' 



^p'{^f + \{Vp'f + u',{p.'^f 



with the changed four-point coupling 

Uo = (det A)"^/2mo : 
with the changed (non-cubic) volume 

d 



y' = []i', = (dctA)-i/V, 



(10) 



(11) 



(12) 



(13) 



(14) 
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L' = Lx: 



^1/2 



(15) 



with a transformed cutoff A' in k' space, k' = A^'^^^Uk, 
and with a critical point at 



(16) 

The temperature variable tq — tqc = ''o ^ ?'oc — '^o^ ''6" 
mains invariant under the transformation ||SJ) and @- 

According to Eq. (|12|l . the bulk critical behavior 
of the anisotropic model H with the coupling uq, Eq. 
((Tj) , can be calculated within the minimally renormalized 
isotropic bulk theory (V' — > oo, A' — > oo) for 7/' with the 
coupling Ug in 2 < d < 4 dimensions |32l |. provided that 
Ug > 0. Specifically, the renormalized quantities of the 
Hamiltonian H' are defined as 



u = II 



-'Z-'zl,A,u', 



r = Z,7^(ro-roJ 
r(3 - d/2)22-'^7r-'^/2(rf - 2)-i and 



with 



(17) 



(18) 



(19) 



(20) 



e = A — d, where S{e) and the Z factors Zi{u' , e) depend 
on £ and u' in the same way as they depend on e and u 
in the standard (A = 1, y — + cx), A — + oo) theory [3^ . 
with an identical fixed point value u'* = u* . This state- 
ment applies also to the field-theoretic functions Cr(u') 
and (^ipi (u') which determine the critical exponents v and 
77. This proves that the critical behavior of H and H' 
belongs to the same universality class in the whole range 
of A where det A > 0. 

Our model, Eq. |71l with A 7^ cgl, can be consid- 
ered as the continuum version of a ip* lattice Hamiltonian 
Hiattice with short-raugc interactions Jij (see, e.g., Eq. 
H5()|l below for a lattice model with a single lattice con- 
stant a). Non-cubic anisotropics may arise either from a 
non-cubic lattice structure or from non-cubic interactions 
on a cubic lattice (as an example see Eq. (|51|l below) or 
from both types of anisotropics. In some range of A 
near A « cqI with cq > 0, Hiatuce and H belong to the 
same universality class. Note, however, that in general 
roc,iattice{Jij]UQ,a) ^ roc(A;Mg, A). 

In order to elucidate the effect of the non-diagonality of 
the anisotropy matrix A we first discuss the bulk order- 
parameter correlation function for T >Tc 



G(x; A,Ug) EE< (^(x)(^(0) >H 



(21) 



where < ... >h means an average with the exponential 
weight . Equations (O, ® and (fT^ imply 



G(x; A, itg) = (det A)-^'^G' {x! ■ u'^) (22) 



where 



G'(x>;,) ^< ^'(x')(/.'(0) >H' 



(23) 



The second- moment bulk correlation length ^'(ug) asso- 
ciated with H' is defined by 



1 rrf'*:E'x'2G'(x>') 

— lini 

2dv'^oo J d'^x'G'(x';u'o) 



1/2 



(24) 



For T > Tc and |x'|/^' < 0(1) the asymptotic scaHng 
form of G'{x.';u'o) reads 4,, 5] for |x'| > A'-^, > A'^^ 



(25) 



with a universal scaling function (f>, a nonuniversal am- 
plitude Ag{uq, A'), and with ^' = Co("o)^ "^j apart from 
corrections to scaling. Eqs. H22|l and H25|l imply 

asymptotically 

G(x; A,^.g) = A'a\ A-^/^Ux f'^'^^d A-^/^Ux | /O 

(26) 

with A'q = ^^(detA) Thus the anisotropy does 

not change the universal structure of the scaling function 
$ but makes the scaling argument of $ and the spatial 
behavior of G anisotropic, even right at T'c(A) (see also 

mi). 

Choosing x = Xiei along the principal direction i,i — 
l,...,d defined by the eigenvector e^, we have (Ux)j 
XiSij and 



G{xie,;A,uo) = A'q{\ \/\ 



l/2N-d+2- 



"$(1 X, I 



where 



e(')(A,^.g)=C^''t 



(27) 



(28) 



are the principal correlation lengths of the anisotropic 
system with the nonuniversal amplitudes 



1/2 f/ 



(29) 



(The amplitudes ^g and ^g'' may depend, in general, also 
on the cutoff.) Their product 



Korr(A) =[|C 
1=1 



(30) 



constitutes the appropriate measure of the correlation 
volume whose shape is ellipsoidal rather than spher- 
ical. This is seen by determining the singular part 
Fs{t]A,uo)/VkBT = fs{t;A,uo) of the bulk free en- 
ergy density f — — limy^oo In J 'Dipe~^ of the 
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anisotropic system. Using Eqs. ifTUI) . ifTT)) . ((T^ and lfTB|) 

we obtain 

/,(t;A,uo) = (detA)-i/2/^(t;<) (31) 

where f'^{t;u'^^) is the singular part of the bulk free en- 
ergy density /' = — limv"-^oo ^ """^ In/ Vip'e'^ associ- 
ated with H', Eq. Together with Eq. 101 for the 
isotropic system, Eqs. (^, and (|2Hl - (|SU lead to 

^hnY fs{t; A, Mo)Vcorr-(A) = (5((i, n) = universal (32) 

which is identical with Eq. JSJ. 

iFiom Eqs. (|26|l - (|29|1 we see that a complete knowl- 
edge of the asymptotic behavior of the correlation func- 
tion G requires the knowledge of the d + I nonuniversal 
amplitudes A'q and ^q'-* and of the d{d — l)/2 nonuni- 
versal parameters characterizing the directions of the 
d eigenvectors e^. For real magnetic materials these 
quantities are unknown as they depend on all micro- 
scopic details. Furthermore, real magnetic materials may 
have lattice structures and anisotropic interactions (e.g., 
ANN, ANNN and third ANN interactions) correspond- 
ing to a nondiagonal matrix A. It is because of the non- 
diagonality of A that both a scale transformation and a 
rotation is necessary and that a simple rescaling of d am- 
plitudes is not sufficient. Clearly two-scale factor univer- 
sality is absent in the bulk correlation functions of such 
anisotropic systems (e.g., metamagnets [s^) although 
they belong to the same universality class as isotropic 
systems (e.g., fluids). 

While the anisotropy does not destroy the univer- 
sality of the scaling function $ of the bulk correlation 
function G (in the non-exponential regime r/^ < 0(1) 
|3|), a fundamental complication arises for confined sys- 
tems since, in general, the principal directions ei of the 
intrinsic anisotropy are totally unrelated to the orien- 
tation of the surfaces of the confining geometry (e.g., 
Li X L2 X ... X Ld rectangular geometry). This intro- 
duces a source of non-universality that cannot be ab- 
sorbed only by a transformation of the lengths Li of the 
confining geometry or of the scaling argument. Within 
our model {Tj), a complete information of this source of 
non-universality requires, at /i = 0, the knowledge of 
d + d{d — l)/2 = d{d -f l)/2 nonuniversal parameters 
(rather than d parameters). Within this model we shall 
show that this implies not only the absence of two-scale 
factor universality but the absence of universality itself 
for all finite-size scaling functions and finite-size ampli- 
tude ratios of anisotropic systems with non-cubic sym- 
metry. In particular, two-scale factor universality and 
isotropy cannot be restored by an anisotropic scale trans- 
formation for confined systems in rectangular geometries 
with a non-diagonal anisotropy matrix A. This is the 
central general result of this paper to be demonstrated 
in the following on the basis of exact results in the large-n 
limit and of one-loop RG results for n = 1, 2, 3. 



First we consider the susceptibility x (psr component) 
of the field-theoretic model ^ above Tc in a finite cube 
with periodic boundary conditions. In the limit n ^ 00 
at fixed uqu it is determined by |34l | 

=ro + 4uoni-'^^(x-'+k.Ak)~' (33) 

k 

with k • Ak = p Aapkakp. The sum runs over k 
vectors with components ka ~ ^ixmaj h^rria = 0, ±1, . . . 
up to some cutoff A. For 2 < d < 4 the asymptotic 
form of the correlation length ^' defined by Eq. (|24|) is 

^ = ^^ii/Ca-'d) with 

= (4u[,nA,ao-V£)'/^'"'^ • (34) 

For large L ^ A^^ and small < i ^ 1 we find the 
asymptotic scaling form for i'/^' ^ ^(1) 

X(i, L- A) = L'^l-g,.^^,{L' li!- A), -ijv = 2 (35) 
with the rescaled length 

i' = L(detA)-i/2'' (36) 
and the normalized anisotropy matrix 

A = A/(det A)!/'^ (37) 
where gcube{x',A.) is determined implicitly by 

- = (4 - d)A-'l,{g^:,^;A) , (38) 



/j(z;A)=/ ds(47r2)-Js^-ip(s,A)e-^''/'*'^', (39) 
"'0 

with 

P(s, A) = (n/sf/^ - J2 e-™-^™" . (40) 

m 

The sum runs over m = (mi, . . . , TOc() with all 

integers TOq = 0,±1,.... For A = 1, gcube{x;l) = 
Qrii.he j.in(x) IS thc kuowu scallug function of the isotropic 
case |3J|. For A 7^ 1, however, gcube{x',A) is nonuniver- 
sal and depends on A in a highly complicated way via the 
inhomogeneous m ^ modes, even after having intro- 
duced the rescaled length L', Eq. H36|l . The effect of these 
modes depends on the orientation of the eigenvectors 
relative to the shape of the confining geometry. In general 
this anisotropy effect cannot be inferred from the knowl- 
edge of finite-size scaling functions of isotropic systems of 
the same universality class and cannot be described by a 
transformation of the argument x of gcube,iso{x) (unlike 
the case for the scaling function $ of the bulk correlation 
function G) or by a rescaling of L. 

Only in the special cases where A = A is diagonal at 
the outset and where the eigenvectors e; happen to be 
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parallel to the edges of the confining cube, the finite- 
size scaling function of the anisotropic system in a cu- 
bic geometry can be reexpressed in terms of the scaling 
function of the isotropic system in a L'l x . . . x L'^ block 

— 1/2 

geometry, L'^ — L\- . Such special cases with a diago- 
nal matrix A are c? = 2 or d = 3 spin models on sc cubic 
lattices with only NN couplings 7^ >/y 0j or 
Jx Jy 7^ Jz ll^, respectively. 

We note that a conclusive answer about the appropri- 
ate way of rescaling the length L cannot be inferred only 
on the basis of the result of x(0, L\ A) at T^, without 
further knowledge. The same statement applies to the 
correlation length A) in a L'^^'^ x 00 cylinder. It 

would always be possible to rewrite x ^-t Tc in the form 



x(0,L;A) = L^/''5, 



the. 



so(0) 



(41) 



with the amplitude 5ciibe,jso(0) of the isotropic system if 
all anisotropy effects are formally absorbed in the length 



L' [.gcufce(0; A)/gc«6e,jso(0)] 



(42) 



But after the calculation of a different physical quantity 
at Tc it becomes obvious that this length L is inappro- 
priate as will be demonstrated in the following. 

Next we present the anisotropy effect on the finite-size 
scaling function of the singular part of the reduced free 
energy density per component in the large-n limit for 
cubic geometry and periodic boundary conditions. For 
L > A-i, < i < 1, L'/e' < 0(1) we find 



/,(t,L;A) =L-'*yen6e(i^7?';A), 



(43) 



In 2 [d-2)Ad 
~Y ^ 2d(4 - d) 



[3(3;; A)] 



-d/2 



8772 



ds 



47r^ 



P(s,A)e-'*/(4^'s) (44) 



where g{x]A.) = gcube(x;A.) is determined implicitly by 
Eqs. (|2Hll - EOJ- The scaling function Ycube{x]A.), in- 
cluding the amplitude l'ctt&e(0; A), is nonuniversal. Only 
on the level of a lowest-mode (k = 0) approximation 
in Eq. (|33|l the explicit dependence of Y^ube on A dis- 
appears. The effect of the m 7^ modes cannot be de- 
scribed simply by a transformation of the scaling variable 
X of Ycube{x] 1) = Ycube,iso{x) of the isotropic case and it 
depends on d{d -I- l)/2 — 1 nonuniversal parameters con- 
tained in A. [Equivalent parameters appear already in 
G, Eq. H26|) .] This holds, of course, also for the rele- 
vant case of general finite n < 00 as can be shown [35j 
within a one-loop RG calculation for the model (|7|) . The 
exact scaling function Ycube{x]A.) for general n remains 
unknown even if the exact scaling function Ycube,iso{x) 
were given for general n and if the exact matrix A were 
given for a special anisotropic system. 



We note that the same rescaled length L', Eq. 136|l . 
is employed in the scaling argument i'/C of Ycube as in 
gcube but not in the leading L^'^ power law of Eq. 
At T — Tc, it would of course be possible to rewrite /s 
in the form 



/s(0,i;A) 



L-^Yr^^ 



cube.isi 



o(0) 



(45) 



with the amplitude Ycube,iso{0) of the isotropic system if 
the anisotropy effect is formally absorbed in the length 



L = L [YcubeiO; A)/YcubeMO)] 



-i/d 



(46) 



This length L differs, however, from the length L, 
Eq. H42() . introduced formally for the susceptibility 
X(0,L;A). 

As seen from our results for x(i, L; A) and fs{t, L; A), 
a possible ambiguity of defining a rescaled length disap- 
pears after calculating the complete temperature depen- 
dence of the finite-size scaling functions of the anisotropic 
system. At the same time such results clarify whether or 
not isotropy can be restored by a scale transformation. 
The exact analytic form of our results (|35|l and H43|l for 
T >Tc unambiguously answers this question for cubic ge- 
ometry and periodic boundary conditions. An extension 
of our results to rectangular LiX L2X ...x Lj, block geom- 
etry fs^l confirms our findings, i.e., even after a rescaling 
of the lengths Li the finite-size scaling functions remain 
nonuniversal for systems with a nondiagonal matrix A. 

We conclude that, for rectangular geometry and pe- 
riodic boundary conditions, finite-size scaling functions 
are, in general, not universally determined only by the 
bulk universality class but do depend on nonuniversal 
parameters in a highly complicated way if the system is 
anisotropic in the sense specified above. In particular, 
within our model ((TJ, if the matrix A is nondiagonal, 
isotropy cannot be restored by a rescaling of lengths 3(| . 

We expect that this conclusion holds also for nonperi- 
odic boundary conditions and for non-rectangular geome- 
tries. For example, we expect that the universality of the 
amplitude u of the "corner" term uL~'^\n.L of the d = 2 
and d = 3 free energy density for free boundary condi- 
tions at Tc 113) nil is not generally valid for anisotropic 
systems. The universality of u was proven in j^^l only for 
isotropic {d — 2) systems whereas in [^^ it was supposed 
to be valid " within a given RG universality class" . Fur- 
thermore, there have been calculations [l3 . 120I [21I 
of edge exponents of anisotropic spin models in wedge 
geometries with free boundary conditions. It was found 
that the anisotropy enters explicitly into the exponents 
and that it was possible to rescale lengths anisotropi- 
cally to bring the problem into an isotropic form. We 
expect, however, that this is, in general, not possible for 
the temperature-dependent finite-size scaling functions of 
lattice systems with edges whose continuum limit yields 
an effective Hamiltonian of the form of Eq. Q with a 
nondiagonal matrix A. 
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We briefly extend our analysis to oo"^"^ x L film ge- 
ometry with periodic boundary conditions, with L being 
the thickness in the d}^ direction. In the large-n limit we 
find 

fsjumit, L; A) = L-"" [(A-i)rfrf)] ""'^ YjumM^) , 

(47) 

where YfUm^iso is the scaling function for the isotropic 
system, with a transformed argument x = i/^'. 



L=[{A-%,)f^L 



(48) 



and where (A"-^)^^ and (A^^)^^ are the d*'' diagonal 
elements of the inverse of A and A, respectively |l32i|. In 
contrast with a nonuniversal amplitude appears at 
bulk Tc and the Casimir amplitude 

A = (d - 1) [(A-l)drf] YfUrnMO) (49) 

is nonuniversal. The simplicity of this anisotropy effect is 
due to the one-loop structure of diagrams contributing to 
the large- ri limit. From finite-size theory at order Ug |40j 
we infer a highly complicated A dependence of fsjum 
for finite n. Furthermore we expect that the amplitudes 
[4ll | and scaling functions '42*1 of density profiles in film 
geometry are nonuniversal for anisotropic systems with 
non-cubic symmetry. More generally, our results suggest 
that the feature of universality in the theory of boundary 
critical phenomena ji^, 143. l45j as well as the notion of 
a "surface universality class" and of " (2 + l)-scalc factor 
universality" need to be reconsidered for the case of 
anisotropic systems. 

It would also be interesting to interpret finite-size 
effects in percolation problems of anisotropic systems 
^3.113 the light of the results of the present paper. 

We illustrate our theory by the example of the Binder 

cumulant ratio [/ for i oo at Tc, Eq. |SJ. We consider 
4 , . . . 

a (f 



Hla 




-si'"' 



(50) 

with an isotropic ferromagnetic interaction Jij = J > 
between nearest neighbors but an anisotropic interaction 
Jij = J' with only 6 (rather than 12) next-nearest neigh- 
bors in the ±(1, 1, 0), ±(1, 0, 1), and ±(0, 1, 1) directions 
on a simple-cubic lattice with a lattice constant a in a 
cube with periodic boundary conditions. It is expected 
that a ferromagnetic critical point exists not only for 
J > 0, J' > but also for J > 0, J' < 0. In the con- 
tinuum limit (a — > 0) this model is reduced to Eq.Q 
with 



(51) 
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FIG. 1: Cumulant ratio 1 — 11(^1) /U{0) vs coupling ratio 
w = J' /{J + 2.7') of the field-theoretic model, Eq. 0, in 
three dimensions for n = 1,2,3 (solid, dotted, dashed lines) 
according to Eqs. 1531 - 15711 . 



where w = J' /{J + 2J') < i and cq ^ 2{J + 2J') > 0. 
(Note that the matrix A is diagonal for a model with 
isotropic NNN interactions.) The positivity of cq requires 
J' > — J/2. For w ^ the eigenvectors are not parallel 
to the cubic axes. The constant cq can be absorbed in 
the bulk amplitude of (|24|l . and does not appear in 



A(m;) = (1 -3^2 + 2^3)- 




(52) 



One of the eigenvalues of A vanishes at Wc — — |, i.e., 
J' = —J/4 (the two other eigenvalues vanish a.t w = 
1, J' = — J). Thus w may vary in the range — ^ < w < ^ 
corresponding to a line of ferromagnetic critical points 
Tc{w) terminating at a Lifshitz point Tciwc) of the (p^ 
continuum model Q (but not necessarily of the iy9^ lattice 
model (|50l) whose line of critical points Tc(w) may end 
at a value of w different from — ^). 

^From a RG treatment of the model ^ within the 
minimal renormalization scheme in three dimensions 
parallel to previous work 46] we obtain U{w) for L oo 
at Tc{w) in one- loop order for n = 1 as 



(53) 



where 



/dss™exp(-iys2-s4) 

^ra(Y) = ^ . (54) 

/dsexp(-irs2-s4) 
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Here the quantity Y depends on w through A.(w), 

Y = -fc |y + A)] + i + 47rfp + /2(P; A)]| , 

(55) 

with 

6 = 144w'*i92(0) , (56) 



f= [247ri/2M'*i/2^2(0)]^/^ (57) 

and u'* = u* = 0.0412 where Ij{z;A) is given by Eq. 
(|39|l . Clearly there is no way of ehminating the comph- 
cated internal dependence on the anisotropy matrix A in 
Eq. H55|l. thus isotropy cannot be restored by means of 
a scale transformation. 

We have also extended this result to general n [35j |. 
While the w dependence is weak for —0.4 < ui < ^ it be- 
comes appreciable upon approaching Wc = — ^ , as shown 
in Fig. 1 for n — 1,2,3. This proves the nonuniversality 
of U{w). Similarly one can derive a w dependence of the 
Casimir amplitude A(w) and of other scaling functions. 
Note, however, that because of the nonuniversal charac- 
ter of U{w) and A(ii;), these quantities may, in principle, 
differ, e.g., for the (d = 3,n = 1) field-theoretic model 
[Eq. lO with the matrix (|5T|l ]. and the (d = 3, n = 1) 
Ising model (with NN and ANNN couplings) even if the 
geometries and the boundary conditions are the same in 
both models. 

This kind of nonuniversal finite-size effect should exist 
near critical points of real systems and should be de- 
tectable in Monte Carlo simulations of d = 2 and d ~ 2> 
spin models. 

Support by NNSFC, by Max-Planck-Gesellschaft, by 
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